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We analyze the turbulence driven particle transport in Texas Helimak (K. W. Gentle and Huang He, Plasma 
Sci. and Technology, 10, 284 (2008)), a toroidal plasma device with one-dimensional equilibrium with magnetic 
curvature and shear. Alterations on the radial electric field, through an external voltage bias, change spectral 
plasma characteristics inducing a dominant frequency for negative bias values and a broad band frequency 
spectrum for positive bias values. For negative biased plasma discharges, the transport is high where the 
waves propagate with phase velocities near the plasma flow velocity, an indication that the transport is 
strongly affected by a wave particle resonant interaction. On the other hand, for positive bias the plasma 
has a reversed shear flow and we observe that the transport is almost zero in the shearless radial region, an 
evidence of a transport barrier in this region. 
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I. INTRODUCTION 

Improvement in magnetically confined plasmas in 
toroidal devices, as tokamaks and stellar ators, is limited 
by the anomalous particle transport at the plasma edge 
[l[ driven by the electrostatic turbulence observed at this 
region (HQ- In the last years, several works have investi- 
gated the turbulence and transport reduction by impos- 
ing an external electric potential that changes the radial 
electric field profile [4|-|6|]. However, more observations 
are required to complete the experimental description 
and identify the proper theoretical concepts that could 
lead to a complete interpretation of the transport con- 
trol achieved by modifying the electric field. Thus, sev- 
eral experiments have been performed in tokamaks and 
stellarators. 

Observations in other type of magnetic confinements 
can also help to expand the experimental description of 
the effects of the radial electric profile on the turbulence 
and the transport. Thus, to improve the experimen- 
tal knowledge of basic turbulence and transport, experi- 
ments have been performed in some devices which confine 
plasmas with selected characteristics of a fusion plasma 
in a simpler geometry and with better diagnostics than 
those possible in major confinement devices @~Ei- 

These basic plasma devices can present a cylindrical 
geometry as the LAPD which was used to study 
the influence of the bias in the turbulence and trans- 
port [ll[ . Other machines with toroidal geometry present 
the magnetic field line curvature in combination with 
plasma gradients as present in fusion experiments. One 
possibility for these plasma toroidal devices is the he- 



limak configuration, present in the BLAAMANN [12|, 
TORPEX tO and Texas Helimak [lj], which has been 
used to study electrostatic instabilities, turbulence, and 
transport (l5l 16j . The helimak configuration is a basic 
plasma toroidal device with characteristics of a fusion 
plasma in a simple geometry 01131111, with a sheared 
cylindrical slab that simplifies the turbulence descrip- 
tion and provides results that can be used to understand 
plasma edge and the scrappe-of layer transport in major 
fusion machines. Altough a stationary equilibrium is not 
expected in this configuration, a turbulent state with 
a peaked radial pressure profile exists |l2| . As the plasma 
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in hclimaks is colder and less dense when compared with 
tokamaks, it is possible to use a large set of diagnostic 
probes. 

The Texas Helimak is also capable of producing states 
of greatly reduced turbulence by biasing and has an in- 
dependent spcctrocopic diagnostic of the plasma flow. 
These features make the Texas Helimak an interesting de- 
vice to study the plasma flow shear influence on the par- 
ticle transport. Some of the analyzed discharges present 
a reversed shear plasma flow and are especially adequate 
to investigate transport barrier onset predicted to this 
kind of flow [13-111 • 

In this article we analyze the turbulence driven parti- 
cle transport in Texas Helimak and investigate how al- 
terations on the radial electric field, through an external 
voltage bias, modify the turbulence and the transport. 
We employ spectral analysis to identify the main changes 
on the power spectra due to the external alterations on 
the radial electric field profile. Moreover, the biasing also 
allows observing the dependence of turbulence spectrum 
and transport on the flow velocity shear. 

Thus, we calculate the transport radial profile and its 
dependence on the electric bias. We observe that the ex- 
ternal bias value changes spectral plasma characteristics, 
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inducing a dominant frequency for negative bias values 
and a broad band frequency spectrum for positive bias 
values. Because of this difference the shots with nega- 
tive and positive bias values are analyzed separately. For 
negative biased shots, the transport reaches a maximum 
where the waves propagate with phase velocities approx- 
imately equal to the plasma flow velocity. We interpret 
this maximum as evidence that the transport is strongly 
affected by a wave particle resonant interaction. On the 
other hand, for positive bias the plasma has a reversed 
shear flow and we observe that the transport is almost 
zero in the shearless radial region indicating an evidence 
of a transport barrier in this region. 



II. EXPERIMENTAL SET UP 

We analyze experiments performed at the Texas He- 
limak 14|, Fig. [T] (a), a basic plasma toroidal device 
located at the University of Texas at Austin. In this 
machine, the combination between the toroidal and the 
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small vertical field creates a helicoidally geometry of its 
magnetic field lines with curvature and shear. The He- 
limak geometry is an approach to a sheared cylindrical 
slab [18| since connection lengths are long enough to ne- 
glect the end effects. Most of these magnetic field lines 
starts and terminate into four sets of four plates located 
at 180° apart at the top and the bottom part of the ma- 
chine. These plates are used as a support to the Langmuir 
probes and to apply external electric potentials (bias) to 
changes the radial electric field profile. In the analyzed 
experiments, the dominant toroidal field is about 0.1 T, 
which addicted with the weaker vertical field create mag- 
netic field lines with « 40m of connection length at the 
middle of the machine (R = lm). The field lines are thus 
helices as shown in Fig. [T](a), spiraling from bottom to 
top. 

Texas Helimak has a vacuum vessel with rectangular 
cross section with external radius, ^external = 1.6 m, 
internal radius, Rmternai = 0.6 m, and height — 2 m. 
For the experiments analyzed in this work, Argon gas 
at 10~ 5 Torr was heated by ECRH with 6kW of power 
inserted by a window located on the inner side of the 
vacuum vessel. The shot duration is up to 20 s and the 
plasma is in a steady state with stationary conditions 
during 10 s, the time interval considered for fluctuation 
analyzes described in this work. 

The diagnostic system count with more than 700 Lang- 
muir probes mounted at the four sets of bias plates. The 
analyzed data were taken by two digitizers, one with 16 
channel and 100 kHz of sample rate and another one with 
64 channels and 7 kHz of sample rate. One of the used 
probe distributions is shown at Fig. Q] (b), where the 
red dots represent the probes used to measure the ion 
saturation current, blue cruxes the probes used to the 
floating potential, and the blue lines indicate the bias 
plate position. This configuration is used to measure the 
turbulence induced radial transport. 
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FIG. 1. (a) Cross section of the Texas Helimak vacuum vessel 
showing the position of the four sets of bias plates (in yellow 
are those used for this work) and a sample of its helicoidally 
magnetic field lines, (b) One of the probe distributions used in 
this work: the red dots represent the probes used to measure 
the saturation current, the blue cruxes the probes used to 
the floating potential and the blue lines mark the bias plate 
position. 



Figure [2] (a) shows the mean density radial profile esti- 
mated from saturation current fluctuations. The density 
profile has a maximum at R « 0.95 m and decreases 
more smoothly on the external side of the peak, which 
corresponds to the low field side of the toroidal magnetic 
field. Figure [2] (b) shows the mean floating potential ra- 
dial profile which presents a minimum at R w 1.1 m and 
a positive gradient in the analyzed radial region, 1.1m 
< R < 1.25m. In this low field side region, the density 
gradient is practically uniform. 

We show in Fig. [3] (a) and (b) examples of evolution 
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FIG. 2. Radial profiles of (a) mean plasma density estimated 
from current saturation and (b) mean floating potential. 



of the normalized saturation current and floating poten- 
tial during a short interval of 0.1s (of the 10 s consid- 
ered for fluctuation analyzes) in a shot without biasing. 
The saturation current signal presented in Fig. |3Ia) was 
measured at R = 1.14m and the floating potential of Fig. 
[3fb) at R = 1.13m (both in the radial interval chosen 
for fluctuation analyses). For the normalization of the 
floating potential presented in Fig. [3]Jb) we consider kTe 
= 16.2 eV, the temperature obtained through the I x V 
electric probe curve. Moreover, the turbulent floating 
potential (red cross) and the ion saturation current (blue 
dot) fluctuations depend on the radial position as shown 
by the radial profile of their turbulence levels shown in 
Fig. 02[c). For each radial position presented in Fig. [3Jc) 
we consider the local average temperature and calculate 
the fluctuation standard deviations (er„ and a v ) and the 
mean density (n) during the whole 10 s of stationary equi- 
librium. To obtain a v we consider the plasma potential 
and the floating potential fluctuations similar. Within 
these approximations, the two curves in Fig. EJc) are 
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FIG. 3. Evolution of normalized fluctuation of (a) ion sat- 
uration current and (b) floating potential in a shot without 
biasing. The saturation current signal was measured at R = 
1.14m and the floating potential at R = 1.13m. (c) Radial 
profile of the density (blue dots) and floating potential (red 
crosses) turbulent level. The lines are guides to the eyes. 

close enough to consider the plasma adiabaticity satis- 



Ill. TURBULENCE CONTROL 

The turbulence level and turbulence driven particle 
transport are controlled by changing the radial profile 
of the electric field through the imposing of an external 
electric potential on some of the 16 available bias plates 
(see Fig. [T](a)). For the shots considered in this paper, 
bias is imposed in four bias plates (two on the top and 
two on the bottom, on both sides of the machine) placed 
in the interval from R=0.86m to R=1.07m, as marked 
in yellow in Fig. [1] (a), near the radial region chosen to 
analyses the transport (l.lm<R<1.25m). 

Figure 2] shows an example of how the external bias 
changes the plasma turbulence. This figure presents the 
frequency spectrogram of the saturation current fluctua- 
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tion for three different bias values. The spectral power, 
showed in grey scale, is calculated by using a windowed 
FFT algorithm and separately normalized for each case. 
The ten seconds of each time series is divided into 200 
time windows with fifty thousand points each. So, for ev- 
ery time window the power spectrum is plotted vertically 
on gray scale. 

For all cases, the equilibrium plasma does not change 
during the whole ten seconds of the analyzed time. The 
spectrum for grounded bias, Fig. |U(a), presents a con- 
centration of power near f 00 Hz with a large width. Com- 
paring Fig. 2](b) obtained for bias = -8V with Fig. |4](a), 
we recognize that the negative bias induces a dominant 
mode at frequency 100 Hz, making the power spectrum 
thinner. On the other hand, in Fig. |U(c) for +8V of bias, 
the spectrum is broader when compared with the unper- 
turbed spectrum. The effect of external bias on the power 
spectra is better observed in Fig. [4] (d) which shows the 
time average of the three windowed power spectra of Fig. 
I (arc). 

The analysis of the saturation current turbulence 
changes induced by the bias, presented in Fig. 01 is ex- 
tended in Fig. [5J where the radial profiles of spectrum 
width (A/) are calculated for several bias values. The A/ 
is estimated by sorting the power spectrum channels in 
crescent order and then computing the number of chan- 
nels that needs to be summed in order to reach half of the 
total spectral power. As shown in Fig. O the external 
bias induces different behavior at the plasma turbulence, 
enlarging the broad band for positive bias and creating 
frequency localized modes for negative bias. Because of 
these different results, turbulence will be separately an- 
alyzed for positive and negative perturbing bias. 

To end the section, we point out that the adiabatic con- 
dition, shown in Figure Etc) for the case without external 
biasing, is still a reasonable approximation for both neg- 
ative and positive imposed bias values. 



IV. RESONANT DRIVEN TRANSPORT 

In this section we use data collected by the probe dis- 
tribution shown in Fig. [TJ (b) to analyze turbulence and 
transport alterations induced by negative bias. 

We consider the spectral analysis introduced in [23[ to 
calculate the particle transport induced by the electro- 
static turbulence. Within this procedure the time aver- 
age transport is given by T —< h ■ Ve >, where n, and 
Ve are, respectively, the density and the E x B electric 
drift velocity fluctuations. As indicated in [23| we take 
h oc I sat and Ve oc E z , , where I sat and E z are, respec- 
tively, the fluctuations of the ion saturation current and 
of the z component of the electric field. 

Applying the introduced spectral analysis, we calcu- 
late the radial transport profiles for several negative bias 
values shown in Fig. [5] (a). In this figure we use the 
same color for positive and negative transport once we 
are interested on the absolute values of the transport. 



This figure shows that the transport depends on the ra- 
dial position and the bias values. Namely, for external 
bias between -5V and -10V transport is high in the region 
around R = 1. 10m, while for bias from -20V to -30V the 
maximum is around R = 1.14m. 

In the analyzed experiments, the bias changes not only 
the estimated particle transport but also the plasma flow 
shear obtained from the spectroscopy radial flow profile. 
The influence of the sheared flow on the transport is pre- 
dicted in (24]|25l ]. A high value of velocity shear is con- 
sidered useful to breaks the macroscopic structures in 
the plasma edge and, consequently, to reduce the trans- 
port. To check this prediction, we present in Fig. [B] (b) 
the radial profile of the velocity shear dV z /dR calculated 
for several external bias values. Comparing Figs. [5] (a) 
and (b), we observe, contrary to the prediction, that the 
transport is high in regions with high shear. Thus, the 
lack of transport reduction may indicates that in the an- 
alyzed region of the Texas Helimak the increase of the 
shear velocity does not reduce the turbulent transport. 

Another effect that may contribute to the radial parti- 
cle transport is the wave-particle resonant interaction, as 
considered in [26l - [28j . Applying this framework to Texas 
Helimak, we introduce a Hamiltonian model and verti- 
cal y and radial x normalized coordinates. In this model, 
the toroidal magnetic field is uniform and the equilibrium 
electric field depends on the radial coordinate x, and two 
coupled drift waves propagate in the vertical direction 
y. Thus, the plasma flow, in vertical direction, is due to 
the E x B equilibrium drift and the chaotic transport, in 
radial direction, is caused by the particle guiding center 
E x B drift. Consequently, the wave driven transport de- 
pends on the sheared plasma flow determined essentially 
by the radial electric field profile. 

To apply the considered model we introduce an almost 
integrable Hamiltonian, with two small wave amplitudes 
such that the system is integrable if only one wave is 
considered: 

H(x, y, t) =$ (x) - Vphix + Aisen(k xl y)cos(k y iy) 

+A 2 sen(k x2 x)cos{ky2[y - {V ph2 - V ph i)]^l) 

One of the results presented in [26[ is that the transport 
is essentially determined by the trapping profile U defined 

as, 



U(x) = J -£-[V E (x)-V ph i] (2) 

In these equations A is the amplitude of each w&ve,V p h 
is the phase velocity, k x is the radial wave number, k y is 
the vertical wave number, Ve is the E x B drift wave 
velocity and Bq is the uniform magnetic field. 

The transport dependence on the trapping profile (U) 
can be understood by analyzing Poincare sections of the 
particle guiding center drift obtained from the considered 
Hamiltonian. 
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FIG. 4. Evolution of saturation current power spectra for three different bias values: (a) ground, (b) -8V and (c) +8V. The 
three mean spectra presented in (d) are taken at the same radial position, R = f .18m. 



In the observed plasma edge turbulence, the U val- 
ues are commonly between and 1. Initially, in Fig. 
we present phase spaces of the integrable system for one 
dominant wave, assuming that U is locally uniform. For 
the case where U is close to zero, Fig. [7] (a), there is 
the formation of periodic islands separated by a scpara- 
trix in form of grid. For this equilibrium space phase, 
the perturbing second wave breaks the separatrix, cre- 
ating chaos and causing a high transport channel along 
the broken separatrix. On the other hand, for a higher 
trapping profile, < U < 1, as in Fig. [3(b), the phase 
space contains resonant islands separated by open lines 
that act as barriers, i. e., impede orbit displacement in 
the radial direction. In this last case, the perturbing sec- 
ond wave breaks the island scparatrices and the barriers, 
creating chaos around the islands. However, the chaotic 
transport is much smaller than the one calculated for U 
near zero. Therefore, we can use the predicted transport 
dependence on the estimated equilibrium values of U to 
explain the transport changes with bias. 

To interpret the observed transport profile alterations 
with bias, showed in Fig. H](a), we present in Fig. [5] the 



radial profile of the velocity difference V e — V p h , which is 
directly proportional to the trapping profile U (x). To 
calculate this velocity difference, we take the electric 
drift velocity Ve as equal to the measured flow veloc- 
ity and estimate the principal wave phase velocity V p h 
by a weighted average phase velocity of the S(k, /) float- 
ing potential spectrum. The comparison between Fig. 
[8] and Fig. [6] (a) indicate that the transport qualita- 
tively changes with U(x) as the bias is varied. Namely, 
in the region where U(x) is close to zero (black regions 
in Fig. 8) the transport is high (red regions in Fig. [S] 
(a)). Thus, our estimations indicate that the transport 
is mainly driven by the wave particle resonance predicted 
in the model considered in this section. 

V. EVIDENCE OF A SHEARLESS TRANSPORT 
BARRIER 

In contrast to shots analyzed in the Section HVl shots 
perturbed by positive bias, considered in this section, 
present broad frequency power spectra without a domi- 
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FIG. 5. Radial profiles of spectral width of fluctuating satu- 
ration current for several positive and negative bias values. 
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FIG. 6. Radial profiles of (a) velocity shear and (b) turbu- 
lence induced radial particle transport for several negative 
bias values. 




FIG. 7. Phase space of the one wave integrable system for 
U=0 (a) and for U=0.5 (b). On this phase space x and y are, 
respectively, the normalized radial and vertical coordinates. 
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FIG. 8. Radial profile of the difference between ExB drift 
velocity and phase velocity for several negative bias values. 
This velocity difference is proportional to the trapping profile 
(IT). 
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FIG. 9. Radial profile of the vertical velocity showing the 
radial position where the shear is zero for bias = +10V 



nant mode. Besides that, for discharges with a positive 
bias, the flow velocity radial profile presents a maximum, 
at R = R*, in the analyzed radial interval. Thus, for 
+10V of external bias, Fig. |9]shows an example of veloc- 
ity radial profile with a maximum at R* = 1.13m. The 
presence of this maximum implies that the plasma has a 
shearless region at R* and reverses its shear flow around 
this radial position. 

Transport studies in reversed shear flows in fluids 
2lll22| and plasmas [28ll20| reveal the presence of a trans- 
port barrier in the shearless radial position that impede 
or reduce the radial transport in the reversed shear re- 
gion. In our work, to examine the effect of the reversed 
shear on the transport during discharges perturbed by 
positive bias, we perform a directly comparison between 
the velocity shear and the electrostatic turbulence in- 
duced transport. 

To verify how the transport depends on the radial posi- 
tion and the applied positive bias values, we perform the 
same spectral analysis introduced in section 4. Thus, we 
show in Fig. [TO] (a) the transport radial profile for several 
positive bias values. In this figure we use the same color 
for positive and negative transport. Analogous to Fig. [6] 
(b) , we show in Fig. [TO] (b) the radial profile of the veloc- 
ity shear dV z /dR calculated for several external positive 
bias values. For all analyzed bias values, Fig. [TO] shows 
that the particle transport is very small where the veloc- 
ity shear is null. This observation can be interpreted as 
an indication of a shearless barrier. 

We can adapt the Hamiltonian model introduced in 
Section HVl to interpret the observed transport reduction 
in the shearless region. Thus, to obtain the Poincare 
section of the guiding center orbits in phase space, neces- 
sary to our interpretation, we consider the nonmonotonic 
trapping profile U (x) determined by the flow profile and 
the average phase velocity. In Fig. [TT] for a bias value of 
+10V, we present one example of such Poincare sections 
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FIG. 10. Radial profiles of (a) velocity shear and (b) tur- 
bulence induced radial particle transport for several positive 
bias values. 

obtained for a U (x) shape (shown in Fig. [TTj) estimated 
from the flow velocity shown in Fig. [5] 

The results illustrated in Fig. [TT1 give an indication of 
a transport barrier located near the shearless position, at 
x « 1.2, where U(x) presents a maximum. Consequently, 
the barrier identified in Fig. [TT]can explain the transport 
reduction in the sherless region (see Fig. [TO] (a)). More- 
over, one must note that, for all positive bias, the two 
white regions in Figs. [TO] (a) and (b) clearly indicate a 
transport reduction in the shearless region. 

VI. CONCLUSIONS 

In this work we analyzed the particle transport in 
Texas Hclimak for external biased discharges with fluctu- 
ations obtained from a large set of Langmuir probes and 
a sheared flow velocity measured by a spectroscopy diag- 
nostic. We showed that the external bias value changes 
spectral plasma characteristics, inducing a frequency lo- 
calized mode in negative bias values and a broad band 
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FIG. 11. Phase space of the integrable system, for a wave 
with the average phase velocity, obtained for a non mono- 
tonic trapping profile U(x). On this phase space x and y are, 
respectively, the normalized radial and vertical coordinates 



spectrum in positive bias values. 

We calculated the turbulence induced particle trans- 
port for several negative bias values and the obtained 
results show that the transport in Texas Helimak is not 
reduced in the high velocity shear regions. On the other 
hand, for these discharges, the transport profile has a 
maximum localized near the region where the Ex B drift 
velocity is approximately equal to the wave phase veloc- 
ity. A particle- wave resonant interaction model is used to 
interpret these maxima transport regions as a resonant 
effect that increases the transport. 

For positive bias values the flow velocity radial profiles 
have a maximum inside the analyzed radial region, so 
the velocity shear is zero where the velocity reaches a 
maximum and reverses its signal near this position. We 
observe that the transport profiles are close to zero in 
the same radial positions where the velocity shear is zero. 
This observation is interpreted as evidence of a shearless 
transport barrier expected in reversed shear flow systems. 

In conclusion, our analysis of the transport in Texas 
Helimak reveal new features that may contribute to im- 
prove the understanding of the particle transport and 
their reduction associated to the barrier onset in plasma 



edge fusion machines. 
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